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Abstract 

We construct a global bifurcation diagram of the plane differential 
system 

x^x{l-x)-xy/{a + x'^), ij = y{5 - fiy / x), 
x{t) > 0, y{t) >0,a>0,S>0,/3>0, 

which describes the predator-prey interaction. 

In [1] was considered the following planar differential systems which describes 
a mathematical predator-prey model 

X = rx{l - -f ) - h^fjy, Q 

Here x{t) > 0, y{t) > — the population densities of the prey and the 
predator respectively and r, s, K,h,m — positive parameters ( for more details 
see [1]). 

After the variables and parameters scaling the system takes the form 

x{t) > 0, y(t) > 0, a > 0, S > 0, f3 > 0. 

System ([2]) was investigated with the help of the qualitative theory of the 
ordinary differential equations in [1] where it was demonstrated some properties 
of the system. For example, it was shown that there is an positively invariant 
region in the phase plane and it was found how many steady states are located 
there and it were found their types and so on. During the investigation the 
authors of [1] studied some bifurcations of codim 1 and 2: the turning point 
bifurcation, the cusp bifurcation, the Andronov-Hopf bifurcation, the Bogdanov- 
Takens bifurcation, the loop bifurcation. A special attention was given to the 
study of periodic solutions. The result of [1] were extended in [2]. 
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In [3] we derived some results which supplemented the investigation from 

mm- 

In this paper we fulfil a complete parametric analysis of system (2) taking 
into account mentioned above and our new results. In particular, we consider 
all bifurcations of system (2) of codim 1, 2 and 3. The result of our study is the 
list of all rough phase portraits of the system under study and we can divide 
the entire parameter space into regions with different rough phase portraits. 
Our results give the base for detailed simulations. 

1. For any parameter values a, (3, S the domain H = {{x , y) E : < x < 
1, y > 0} is w-invariant set of The boundary of the domain 11 contains two 
steady states of points (0,0) and (1,0). 

For all parameter values the rest point (1, 0) is a hyperbolic saddle with its 
outgoing separatrix passing into 11. 

The origin (0, 0) is a non-analytic steady state, we change the time t — tx. 
After that ([2]) takes the equivalent form 

X ^ x\l - x) - ^ = P2{x,y) + . . . , 

y = Sxy - Py"^ = Q2{x,y), I® 

x{t) > 0, y{t) > 0, a > 0, (5 > 0, /3 > 0. 

We denote the derivation on r with the point again. P2{x,y),Q2{x,y) are 
homogeneous polynomials of degree two, dots replace terms of degree more than 
two. 

For ([3|) the origin is a complex analytic steady state, the matrix of the linear 
part contains zero elements only. 

As it is proved in ^4J in our case we have the following 

1) if (5 < 1 the neighborhoods of the origin inside R^ contains a unique 
hyperbolic sector with the coordinate axes as its boundary; 

2) if (5 > 1 the neighborhoods of the origin inside i?^ contains one hyperbolic 
sector and one parabolic sector with a common boundary which is the trajectory 
outgoing from the origin with the angle § = atan^^. The second boundary 
of the parabolic sector is the x-axis. 0-trajectories have the asymptotic y = 
kx\k > 0. 

Phase portraits of system (2) near the origin and the infinity are shown in 
the Fig. 1. 

2. To fulfil the complete parametric analysis we first need to study all 
codimension one bifurcations in system (2) to which correspond hypersurfaces 
in the parametric space dividing the parametri space into regions within which 
the system has distinct rough phase portraits. For parameter values on this 
surfaces the system has one from the following singularities 

I. The turning point (the saddle- node). 

II. The fine focus of the first degree. 

III. The double limit cycle. 

IV. The seaparatrix loop with the non-zero saddle value. 

The local bifurcations of codim 1 for I, II can be described by explicit for- 
mulas 
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Fig. 1. 



Let {xo,yo) is a steady state of (2). Then the fohowing algebraic system is 
vahd 

xl-xl + {a + S/P)xa-a = 0, yo = {S/P)xo, 

which has from one to three solutions inside our domain 11. 

Let J is a matrix of a linear part near the steady state. We have 



F 



Xq 



- S/p)xo 



0, yo 



detJ = fci(-a^/3 + 20^/3x0 + 2aSxo - 2a(3xl + AajSxf^ 



tr J = k2{aS + axo 



2a;g 



■Sxl 



(3xt + 2Pxl) 



54)^ ki > 0, 



fc, >0. 



L After the elimination of x conditions F = 0, det J = isolate the bifurca- 
tion set T in the parameter space. The set corresponds to the double steady 
state. T is a piecewise smooth surface and it may be defined as follows 



2t\ P 



2t(1 - r) 



The surface is sketched in the Fig. 2. The points inside the angle correspond 
to the existence of three nonhyperbolic steady states in system (2) (one saddle 
and two antisaddles) and the points outside the angle correspond to the existence 
of a unique steady state (the node or the focus). The points at the surface T 
(but not at TT) correspond to the existence of one double saddle-node and one 
node or focus. The points at the curve TT correspond to the existence of one 
steady state of the multiplicity 3. 

Here and below we examinated the corresponding non-degenerate conditions 
of the bifurcations with the help of results from [2] and our results. 

IL After the elimination of xq conditions F = 0, tr J = 0, det J > isolate 
the bifurcation set H corresponding to the existence of the fine focus in system 
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Fig. 2. 



(2) (the eigenvalues of the matrix of a hnear part are pure imaginary). This 
bifurcation set is a smooth surface with the boundary (detJ = 0). It may be 
defined as foUows 

At aU cases the steady state (xojyo) is given by the expressions 

5 

xo =T, yo = -T. 

Surface H has a comphcate structure. Therefore we do not sketch it herein 
In Fig. 3 it is shown schematicaUy sections of surfaces iJ, T with the hyper- 
planes 5 = c for various values of c. On this sections the surfaces are shown as 
the curves which we label as H, T again. About other notations at the picture 
we explain below. 

Curve T has a cusp. Dy dotted line it is shown the section of an extension 
of surface H for which F = 0,trJ = 0, det J < 0. For corresponding parameter 
values there is a saddle with zero sum of eigenvalues of matrix J. 

In the Fig. 4 it is shown the section of surfaces T by the hyperplane 
5 = .25. We see that curve H has the selfintersection point which is in the 
domain where system (2) has three steady states. 

The conditions about double cycles (III) and separatrix loops (IV) can not be 
describe by explicit formulas in our case. To construct surfaces of this nonlocal 
bifurcations we consider local bifurcations of more high codimensions (two and 
three). Such a local bifurcations are listed in the table. 



^ In [7] there is a picture with Andronov-Hopf bifurcation surface which is analogous to our 
surface H. 
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Table 1: The list of local bifurcations 



Bifurcation 


Notation 


codim 1 




The turning point 


T 


Andronov-Hopf bifurcation 


H 


codim 2 




The cusp 


TT 


Bogdanov-Takens bifurcation 


HT 


The generalized Hopf bifurcation 


GH 


codim 3 




The focus singularity with nilpotent linear part 


HTT 


The cusp singularity with nilpotent linear part 


DHT 



The listed local bifurcations may be described by explicit analitic expres- 
sions. 

If there is a steady state of multiplicity 3 in (2) then we have a solution of 

equation F = with the same multiplicity. In the parameter space correspond- 
ing set is curve TT at the surface T which is given as follows 

a = 1/27, (3 = 275/8, <5 > 0. 

The bifurcation set HT is given by the condition Ai = A2 = (F = det J = 
trJ = 0). The set is a smooth curve and may be described as follows 
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Fig. 4. 



The bifurcation set GH is given by the condition F = tr J = = 0, det J > 
0, where li is the first Lyapunov value of the fine focus. The non-degenerate 
condition of this codim 2 bifurcation is the condition ^2 7^ where I2 is the 
second Lyapunov value. 

With the help of the algorithm of Poincare [9l [10] we derived symbolic ex- 
pressions for the first and second Lyapunov values of the system 



These expressions are 



ao2aii+ai2 -011020 + 3030- 2002602 + 3&o3 + ^02^ii + 2a2o620 + &ii&20 + &2i 



6 



h = 10ao2^aii + 30ao2ao3aii — 9ao4aii — 10ao2^ai2 — 15ao3ai2 — 15ao2ai3+ 
9ai4 + 97ao2^aiia2o + 27ao3aiia20 — 72ao2ai202o — 21ai3a2o + 176ao2aiia20^ 
-80ai2a2o^ + 89aiia2o'' + 15ao2aiia2i - 12oi2a2i + 18aiia2oa2i - 3aiia22- 
75ao2^a3o - 18ao3a3o - 6aii2a3o - 288002020030 - 267a2o^a3o - 27a2ia3o- 
9002031 - 27a2oa3i + 9032 + 15aiia4o + 45a5o + 20ao2^6o2 + 60002003^02- 
18004602 + 34oo2aii^6o2 - 28011012602 + 174ao2^02o6o2 + 24003O20602+ 
31oii^02o6o2 + 208002020^602 + 18020^602 - 3O020O21602 + I2022602- 
63oii03o6o2 + 66040602 + 121oo20ii5o2^ - 71012602^ + 68011020^02^- 

174030602^ + 106002602^ - 18020&02^ - 30oo2^6o3 - 45 003^03 - 9oii^5o3- 

246ao202o6o3 - 282o2o^feo3 - I8021603 - 75oii6o2&o3 - 1596o2^foo3 - 6O002604- 
66020^04 + 45^05 + 31oo2^aii6ii + 6003011611 - 3I002O12611 - 6013611 + 
IIO002O11020611 - 7O012020611 + 79oiiO20^6ii - I2O002030611 - 237020030611 + 
52002^602611 - 3003602611 - 3011^602611 + 56002020602611 - 104020^602611- 
6021602611 - 25oii6o2^6ii - 53602^611 - IO8002603611 - 225020603611- 

2 2 2 2 2 

15604611 + 2I002011611 - 2I012611 + 2I011020611 - 63o3o6ii -4002602611- 

85020602611^ - 60603611^ - 18602611^ + 15oo20ii6i2 - I5012612 + 12011020612- 

36030612 - 3O020602612 - 27603612 + 27602613 + 25002011^620 - 25011012620- 
50002^020620 - I2003O20620 + 22011^020620 - 192002020^620 - 178020^620- 
6002O21620 - 36020O21620 + 6022620 - 54011O30620 + 6O040620 + IIO002O11602620- 
66012602620 + I6011020602620 - I38030602620 + 120002602^620 - 136020602^620- 
66011603620 - I8O602603620 - 25002^611620 - 6003611620 - 3011^611620- 
I82002020611620 - 265020^611620 - 3021611620 - 38011602611620 - 104602^611620- 
43002611^620 - 124020611^620 - 18611^620 - 6002612620 - 36020612620 - 3611612620+ 

2 2 2 2 

9613620 + 25002011620" - 25012620 - 16011020620 - 30030620 + 50002602620 - 
102020602620^ - 75603620^ - 13011611620^ - 61602611620^ - 20020620^ - 10611620^- 
25002^621 - 6003621 - 3011^621 - 102002020621 - 107020^621 - 3021621- 
16011602621 — 44602^621 — 43002611621 — 82020611621 — 18611^621 — 6612621— 
13011620621 - 36602620621 - 10620^621 - 6002622 - 12020622 + 3611622 + 9623+ 
12002011630 - 12012630 + 15011020630 - 9030630 + 24002602630 + 24020602630- 
36603630 + 6011611630 + 9602611630 + 24020620630 + 12611620630 - 3621630+ 
6011631 + 21602631 + 15620631 + I8020640 + 9611640 + 9641. 

Let us note that our expressions coinside (with an accuracy to a positive 
factor) with expressions derived by other authors (see ^ for example). 

In our case we have for system (2) (with an accuracy to a positive factor 
again) 

Zi = 6S-llS^+8S^-2S'^+4:Xo-396xa+^26'^xo-US^xo-24:xl+ 
+725xl-33S^xl++S6xl-335xl-12xl. 

In Fig. 5 it is shown the curve ^i = and the hne xq + S = 1 (the last 
one describes the boundary of surface H). We see that bifurcation set GH is 
a smooth two-component curve. The second Lyapunov value is negative in all 
points from GH (we tested this fact by numerical calculations). Therefor the 
complex focus in system (2) has its multiplicity no more than two and the focus 
is stable if its multiplicity is equal to two. 
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Bifurcations of codim 3 may be detected checking the non-degenerate con- 
ditions for bifurcation of codim 2. 

In [5] it is proved that for parameters 

a = 1/27, (3 = 9/4, S = 2/3 

in the system there is a bifurcation "the focus singularity with nilpotent Unear 
part" ( HTT). 

It is proved also in this paper that for 

a = -45 + 26\/3, /? = (5 + 3\/3)/4, 6 = Vs - 1 

in the sysyem there is a bifurcation "the cusp singularity with nilpotent linear 
part" ( DTT). 

Parametric neighborhoods of listed local bifurcations are well Known and we 
know about sets of non-local bifurcations there. For example, with bifurcations 
HT it takes place non-local separatrix loop bifurcation L\ (Lr), with bifurcation 
HTT it takes place non-local double cycle bifurcation C and so on. 

So we have information about the existence and construction of the non- 
local bifurcation surfaces in regions of the parameter space near sets of local 
bifurcations. More detailed study of the surfaces may be fulfilled with the help 
of the qualitative theory of differential equations and numerical tests. 

As a result of our investigation we obtain a global bifurcation diagram. 

We do not describe here the diagram and we do not describe the process of 
its construction because this description repeats analogous investigations from 
other papers (see [71 [TTl [T2] and references therein) . 

To give an idea of the construction of our bifurcation diagram and to summa- 
rize our results we give the "most complex" section of the bifurcation diagram 
by the plane 6 — c for value of c near 0) (Fig. 6). 
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In the section the parameter space is R\ with the coordinates a, j3, codim 
one bifurcation sets are curves and codim two bifurcation sets are points with 
the same notations as corresponding bifurcations. 

In our figure we sketch only disposition of bifurcation curves and points 
without foUowing the scale because they are very close to each other sometimes 
and it is not possible to draw their full-size images. 

In Fig. 7 are given schemes of rough phase portraits of system ([2|l. 

We show only trajectories within the set 11 and we do not sketch the bound- 
ary steady states. Stable steady states are depicted by black circles, ubstable 
steady states are depicted by white circles, saddles are depicted by half-black 
circles. Unstable limit cycles is depicted by dashed curves. 

We do not show the phase portrait for regions 7'-14' of the bifurcation 
diagram. These portraits may be obtained from the portraits for regions 7-14 
by interchanging two steady state which are non-saddles. 

To complete the list of all phase portraits we show schemes of trajectories 
for two regions 15-16 of the global bifurcation diagram which are absent in the 
section. 

We are thankful to professor Dongmei Xiao (Shanghai Jiao Tong University) 
for discussions and remarks. 
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